Substrate coherency-driven octahedral rotations in perovskite oxide films 
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We perform first-principles density functional calculations to explore the role of substrate proximity 
effects on the octahedral rotation patterns in perovskite oxide superlattices. With cubic perovskite 
SrFeOs as our model film and tetragonal SrTiOa as the substrate, we show that in most cases the 
substrate octahedral rotation patterns propagate into the film across the heterointerface. We also 
identify elastic boundary conditions for which the enforced structural coherence induces atomic 
displacement patterns that are not found in the bulk phase diagram of either individual constituent. 
We suggest that such substrate coherency-induced octahedral texturing of thin film oxides is a 
promising approach for tuning the electronic structure of functional oxide thin films. 



The use of substrate-induced bi-axial strain to modify 
the properties of epitaxial thin films has been demon- 
strated for a wide range of phenomena and materials, 
including mobility in semiconductors, Curie temperatures 
in ferromagnets and ferroelectric polarizations in complex 
oxides [1]. While the main effect of strain in strongly 
covalently bonded materials such as semiconductors is 
to change the bond lengths, the flexible corner-sharing 
networks of oxygen polyhedra in complex oxides provide 
additional routes for accommodating substrate-induced 
changes in lattice parameters. Established mechanisms 
include modifications of the oxygen polyhedral tilt pat- 
terns, and the formation of twin (change in orientation of 
long and short axes) or antiphase (variation in phase of 
polyhedral rotations) domains. 

Recently it is been suggested that, in addition to chang- 
ing the lattice parameter of a complex oxide film, the 
presence of a heterointerface could alter the relative sta- 
bility of polyhedral tilting patterns in both the film and 
substrate through proximity effects |2H5]. While theoreti- 
cal studies have shown that strain-induced competition 
between polyhedral rotation modes and other lattice dis- 
tortions is crucial in determining the functional properties 
of complex oxides |5HTU] , very little is known about the 
extent to which substrate proximity modifies polyhedral 
tilt patterns. This is in part due to difficulties in obtain- 
ing high precision measurements of oxygen positions in 
superlattice and thin film interfaces [TTJ [T^] . In this Let- 
ter, we use density functional theory (DFT) to calculate 
explicitly how the structural distortions of a substrate 
affect the atomic structure and properties of a coherent 
film. 

We take perovskite-structured SrFeOa/SrTiOa as our 
model system, chosen for its continuous A-site sublattice, 
absence of polar discontinuity, and simple oxygen octa- 
hedral tilt patterns: SrFeOs has the ideal cubic Pm3m 
perovskite structure down to the lowest temperature stud- 
ied (4 K) [13], and the ground state I4/mcm phase of 
SrTiOa (which is a widely used substrate) has a single 
octahedral instability with respect to the cubic phase 
[T4] that condenses below ^105 K. First, we investi- 
gate the effect of heterostructure periodicity in symmet- 



ric (SrTi03)„/(SrFe03)„, n — 1 . . . 5, and asymmetric 
(SrTi03)„/(SrFe03)m, n=1...3,m = 1...3 superlat- 
tices. We find that the octahedral rotations from the 
SrTiOs substrate propagate into the first two interfacial 
SrFeOs layers; for highly confined ferrate layers additional 
electronically-driven lattice instabilities occur that are 
not observed in bulk SrFeOa. Then we show that these 
octahedral and electronic lattice instabilities cannot be 
induced in SrFeOa using bi-axial strain alone, indicating 
that substrate coherency and confinement play a critical 
role in determining the heterointerface atomic structure. 

Our DFT calculations are performed within the local 
spin-density approximation + Hubbard U (LSDA-I-/7) 
method using the Vienna Ab initio Simulation Package 
(VASP) [m [TB], with the Dudarev method [T^ for the 
Hubbard correction to the LSDA exchange-correlation 
functional. We use an effective U parameter of 6 eV on 
the Fe d orbitals, and impose ferromagnetic spin order 
in the SrFeOa; this approach is consistent with earlier 
first-principles calculations [THj . Complete details of our 
calculations are reported elsewhere [T^ . 

We first calculate how freezing in the antiferrodis- 
tortive (AFD) a°a'^c~ octahedral tilt pattern of the low- 
temperature phase of SrTiOs [^Hj affects the total energies 
of the two bulk materials which comprise the superlat- 
tices. Figure flla) shows the calculated total energies 
as a function of increasing amplitude of octahedral rota- 
tion angle 9, referenced to the un-rotated states at the 
L(S)DA equilibrium lattice parameters. As expected for 
bulk SrFeOa, we find that the cubic perovskite structure 
is stable with respect to the octahedral rotation mode. 
In contrast, the oPaPc~ tilt pattern is energy lowering for 
tetragonal SrTiOs, consistent with the low-temperature 
experimental structure |30j. 

We next explore whether heterostructuring SrFeOa lay- 
ers with SrTiOs in different period superlattices induces 
AFD FeOg instabilities in the ferrate layers. We exam- 
ine symmetric and asymmetric superlattices, by stacking 
(n, m) five-atom perovskite cells along the c-axis, with 
the in-plane periodicity increased to be commensurate 
with the AFD rotations. We keep the in-plane lattice 
constant fixed to the LDA equilibrium value of 3.86 A to 
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FIG. 1: (Color) Energy versus rotation angle B of the cP<^c octahedral tilt for cubic SrFeOs and tetragonal SrTiOa (a). The 
same mode homogeneously frozen into (b) asymmetric (SrTi03)n/(SrFe03)m and (c) symmetric superlattices (right). 



simulate epitaxial growth on tetragonal SrTiOa (SrFeOa 
is under a L(S)DA theoretical +1.75% bi-axial strain). 
The out-of-plane lattice constant of each superlattice is 
then fully relaxed. 

We first introduce the same octahedral rotation pattern 
as before uniformly into each layer of superlattice. The 
evolution in the total energy for the asymmetric (n, m) 
and symmetric (n, n) superlattices with increasing mag- 
nitude of Q is shown in Figure flTb) and (c). For the 
asymmetric superlattices we find that when SrTiOs com- 
prises the majority of the superlattice, a homogeneous 
rotation of all octahedra is favored over the unrotated 
configuration. In contrast, the symmetric superlattices, 
and those with a greater fraction of SrFeOa layers are 
weakly, if at all, susceptible to the same rotational modes. 

Since the homogeneous rotation throughout the su- 
perlattice is artificial, we next remove the homogeneity 
constraint and allow full relaxation of the oxygen posi- 
tions in each layer within the symmetry of the (Pc^^c" tilt 
pattern. We initialize the oxygen positions in each super- 
lattice to configurations corresponding to the energy well 
minima shown in Figure fl] (For homogeneous rotations 
that were found to be stable, we start with Q = 2°.) The 
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FIG. 2: (Color) The layer-by-layer resolved octahedral ro- 
tation angles (S) for the asymmetric (a) and symmetric (b) 
superlattices. The magnitude of the SrTiOs AFD octahedral 
rotations rapidly decreases into the SrFeOs interfacial layers. 



results of the structural relaxation are shown in Figure 
[2| where we plot the layer-by-layer resolved rotations d 
about the c axis for the different period superlattices. We 
find that in all cases, the octahedral rotations remain 
antiferrodistortive and that the SrTiOa substrate exhibits 
the a^t^c' tilt pattern. Interestingly, AFD octahedral 
rotations are induced in the interfacial SrFeOa layers, 
with the magnitude of the octahedral rotations decaying 
exponentially from the interface into the center of the 
SrFeOa slab. The layers in the center of each slab are 
close to their respective bulk calculated values, indicated 
by the broken lines, although the highly asymmetric (5,1) 
superlattice shows a 9.3% enhancement in Q at the center 
of the SrTiOg slab. 

We next fully relax the atomic structure of each super- 
lattice by removing the symmetry constraint imposed by 
the (^(^c~ octahedral tilt pattern. The asymmetric su- 
perlattices in which SrFeOa is the majority component do 
not show any considerable changes in the atomic displace- 
ments: there is only a small decrease in the magnitude of 
the (^(£'c~ octahedral tilt at the interface due to small 
changes in the apical cation-oxygen bond lengths. We also 
found similar minor changes in the atomic structure for 
the symmetric superlattices (n, n) with n > 2. Therefore 
our earlier conclusion that only the first two ferrate inter- 
face unit cells are modified by the octahedral rotations 
found in the substrate remains for these cases. 

Drastically different behavior is found, however, for 
superlattices with single unit cells of SrFeOa: in the (1, 1), 
(3, 1) and (5, 1) superlattices we find that, in addition to 
the original octahedral rotation pattern, the FeOe octahe- 
dra also exhibit Fe-0 bond length distortions associated 
with electronic instabilities. In each case, either a Jahn- 
Teller (JT) distortion, which produces two short and two 
long equatorial Fe-0 bonds, or a "breathing" distortion 
that makes a uniform dilation and contraction of the 
Fe-O bonds is found to be stable. Note that, for the 
single ferrate layer heterostructures, we obtain JT- and 
breathing-distorted FeOe octahedra even when we dis- 
able the oxygen octahedral rotations in our calculations. 
This indicates that these electronic instabilities are the 
consequence of the quantum confinement of the ferrate 
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FIG. 3: The effect of bi-axial strain on the mode stiffness 
of the rotational patterns in homogeneously strained SrFeOa 
films. For all cases, [001] orientated substrates are considered, 
and the out-of-plane lattice constant is chosen to conserve 
the experimental SrFeOa volume after fixing the in-plane 
lattice parameters to match the following substrates: SrTiOa 
(STO), LSAT= (LaA103)o.3(Sr2AlTa06)o.7, LaAlOs (LAO), 
and YaAlOa (YAO). 



layer, which increases the susceptibility to electron local- 
ization and in turn enhances orbital degeneracy-lifting 
instabilities [11] . In the unrotated case, the energies of the 
JT- and breathing-distorted structures are energetically 
equivalent within the resolution of our DFT calculations. 

While the JT- and breathing distortions are not induced 
by the octahedral rotations, when both rotations and JT- 
or breathing distortions are allowed, we are able to resolve 
distinct ground states, in which the octahedral rotations 
cooperate with a specific electronic instability. In the 
(3,1) superlattice, for example, the ground state consists 
of a large Jahn- Teller bond length distortion of 0.05 A, 
in combination with a completely different octahedral 
rotation pattern - the well-known perovskite a^a^c^ 
GdFeOa tilt pattern - that does not occur in either of 
the parent compounds. If we instead enforce the usual 
a'^a°c~ tilt pattern of SrTiOs, the superlattice exhibits a 
breathing distortion with Fe-0 bond length differences of 
6.6% between the two inequivalent FeOe octahedral sites. 

Finally, we check whether our finding of octahedral ro- 
tations within the SrFeOa interfacial layers can be repro- 
duced using strain alone, or whether it requires substrate 
coherency. To explore the interaction between strain and 
the octahedral rotations on the SrFeOs layers we remove 
the SrTiOa substrate from our calculations and simu- 
late homoepitaxially strained SrFeOs, with the bi-axial 
constraint imposed by enforcing equal in-plane lattice 
parameters. In addition to the ground state a^a°c~ tilt 
pattern, we also examine the strain dependence of the 
in-phase a'^a'^c^ rotation about the c-axis, the a~a'^c^ 
AFD tilt pattern about the a-axis that lies in the epitax- 
ial plane, and the Jahn- Teller and breathing distortions. 
Note that similar studies were performed previously for 
SrTiOa [H EH [22] ; we do not repeat those calculations 
here since our SrTiOa substrates are strain free. 

In Figure |3] we plot the calculated mode frequencies of 
the octahedral rotational patterns for a range of strain val- 



ues corresponding to typical substrates. Real frequencies 
indicate that the cubic lattice is stable and that the mode 
does not spontaneously condense. Importantly, we find 
that all the frequencies are real for strain corresponding to 
coherency on SrTiOs, indicating that the rotations in our 
SrFeOs/SrTiOa heterostructures require the actual pres- 
ence of the heterointerface, not only its associated strain. 
Under > 1% compressive strain, we find an unstable AFD 
rotation which has been reported previously for the case 
of LaAlOa |23] and provides a low energy route to reduc- 
ing the in-plane lattice parameters without significantly 
shortening the bond lengths. Interestingly, the aPa'^c'^ 
tilt is less sensitive to strain and is likely due to the rect- 
angular planar coordination of the Sr-site, which leads 
to less energy stabilization from strain-induced covalency 
modifications over the distorted tetragonal coordination 
available in the a°a°c- tih pattern [21]. The a-a°c° AFD 
tilt also softens with increasing tensile strain to accom- 
modate the reduction in the out-of-plane lattice constant 
as the in-plane lattice parameters are elongated. [The 
Jahn- Teller and breathing distortions (data not shown) 
remain metallic and are disfavored for all strain states.] 

Finally, since we define strain relative to the LSDA 
equilibrium volumes in our calculations, we examine the 
dependence of the SrFeOs lattice instabilities on unit cell 
volume. We find that all phonon frequencies are real 
(indicating no instabilities) for lattice constants within 
±2% of the experimental value (oq = 3.851 A). Also, 
the choice of U does not significantly alter the phonon 
dispersions. 

We now discuss how appropriate choice of the antifer- 
rodistortive octahedral rotations in the substrate may be 
used to select particular latent lattice instabilities in a 
functional thin film. As we showed earlier, if the octa- 
hedral tilt pattern in the film can be modified, it may 
be possible to switch between the types of electronically- 
driven structural instabilities that couple to the rotational 
modes. This is shown schematically in Figure |4] for the 
(SrTi03)„/(SrFe03)i superlattices. When the superlat- 
tice exhibits the orthorhombic Pnma octahedral tilt pat- 
tern a~a~c^^ the Jahn- Teller distortion is activated due 
to the orbital degeneracy in the highly confined SrFeOa 
layers. On the other hand, the a'^a'^c~ tilt pattern is 
found to occur with the octahedral breathing distortion. 
The experimental situation likely depends on the different 
temperature scales associated with the octahedral and 
electronic-driven structural distortions; since the rotations 
condense at higher temperatures, particular substrates 
may be chosen to be compatible with certain symmetry 
electronic-instabilities that condense at lower tempera- 
tures. In the case of these ferrate/titanate superlattices, if 
a substrate that shows orthorhombic or tetragonal octahe- 
dral rotations is used, it may be possible to strictly select 
orbitally- or charge-ordered electronic configurations in 
the ferrate layers by heteroepitaxy. 

An interesting question is whether such strong coupling 
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FIG. 4: Illustration of how the same high symmetry reference 
superlattice can support either the electronically-driven Jahn- 
Teller or breathing distortions by the condensation of different 
symmetry octahedral modes with temperature. The choice 
of a substrate based on the octahedral rotational properties 
may substitute for temperature to control such phases under 
standard temperatures. 



between the octahedral tilt modes and electronic instabil- 
ities is present in bulk perovskite oxides. Such coupling 
is active in the isoelectronic d mangantite compound 
LaMnOa, where the orthorhombic a^a^c^ tilt pattern 
combines with a long-range ordering of Jahn- Teller dis- 
tortions, in preference to the breathing distortion, to 
induce a metal-insulator phase transition [25]. The oc- 
tahedral rotations and La displacements are critical to 
stabilizing the ground state electronic configuration and 
the metal-insulator transition pSll27j . In fact, without 
the tilting the t2g and eg states are forbidden to mix by 
symmetry and the orbitally ordered ground state imposed 
by the Jahn- Teller distortion does not occur [28J. Since 
orbitally-degenerate perovskite oxides support strong cou- 
pling between lattice modes of different origins, growth 
on substrates with different octahedral instabiUties offers 
a new route to tune the electronic properties of these 
materials. 

In summary we used density functional calculations 
and the model SrFeOs/SrTiOa superlattice system to 
show that thin film perovskite heterointerfaces are not 
only sensitive to the elastic strain at the interface, but 
also to the octahedral rotation patterns present in the 
substrate. We found that the a''^a'^c~ tilt pattern in the 
SrTiOa substrate propagates into the SrFeOa film inter- 
face layers. Interestingly, we find that the substrate not 
only affects the rotation patterns of the film, but that the 
entire heterostructure can adopt a tilt pattern that is dif- 
ferent from that of the parent substrate or film structure; 
this may partly explain the recent anomalous structural 
phase transitions observed in SrTiOa substrates used to 
support manganite/cuprate superlattices [4]. We point 
out that the SrTiOs a^a^c~ tilt pattern studied in this 
work is somewhat weakly correlated along the c-direction 
and is the reason for the short penetration depth across 
the interface. Stronger coupling of substrate coherency- 



induced octahedral rotations can be obtained by choosing 
substrates which show in-plane rather than out-of-plane 
octahedral rotations patterns with orthorhombic or mono- 
clinic space groups. Alternatively, control of the substrate 
vicinal cut may be used to control the types of octahedral 
rotational antiphase domain boundaries induced by the 
substrate. The ability to control octahedral rotation tex- 
tures present in thin film via substrate proximity effects 
affords an additional parameter for tuning the atomic and 
electronic structure of functional oxide films. 
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